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I. INTRODUCTION 



Effective field theory calculation of second post-Newtonian binary dynamics 

James B. GilmorcQ and Andreas Ros!0 
Department of Physics, Yale University, New Haven, Connecticut 06520, USA 
(Dated: February 10, 2009) 

We use the effective field theory for gravitational bound states, proposed by Goldberger and 
Rothstein, to compute the interaction Lagrangian of a binary system at the second post-Newtonian 
order. Throughout the calculation, we use a metric parametrization based on a temporal Kaluza- 
Klein decomposition and test the claim by Kol and Smolkin that this parametrization provides 
important calculational advantages. We demonstrate how to use the effective field theory method 
efficiently in precision calculations, and we reproduce known results for the second post-Newtonian 
. . . order equations of motion in harmonic gauge in a straightforward manner. 
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In the last two decades, significant progress has been made towards the detection of gravitational waves (GWs) 
via laser mtcrferometry. Currently, the ground-based experiments LIGO [l|, VIRGO [1], GEO [1], and TAMA [j] 
are actively searching for GWs Moreover, the proposed LISA experiment 0, due to be the first space-based GW 
detector, will search for GWs in a complementary frequency band to the ground-based experiments and is expected 
to achieve high event rates at an unprecedented signal-to-noise ratio Q ■ 

A particularly interesting source of GWs, which is expected to be detected, is the compact binary system undergoing 
bJQ, coalescence, with neutron star (NS) and/or black hole (BH) constituents. Current experiments have yet to detect 
the binary inspiral signal. However, Advanced LIGO an upgrade of LIGO scheduled to come online in 2014, may 
£NJ ■ allow for routine detection of such events. This is due to a ~ 10-fold increase in sensitivity over LIGO, which will 
in turn result in an increase of the accessible event rate by a factor ~ 1000. Current estimates for the number of 
expected NS/NS, BH/BH, and BH/NS events in Advanced LIGO are roughly 10 - 100, 1 - 500, and 1 - 30 per year, 
fT} ■ respectively HQ. 

All three stages of the binary coalescence, inspiral, merger, and ringdown, are potentially detectable. The inspiral 
phase, where the characteristic orbital velocity is v 2 -C 1 (in units where c = 1), can be computed analytically using 
an expansion in v 2 ~ Gm/r. The merger is computed numerically llj, and there has been significant recent progress 



in this area 12 1. The ringdown can be treated analytically using quasinormal modes 

The perturbative calculation of the inspiral phase has been performed with a variety of methods Because 
of the phase evolution of the inspiral signal and the ability to measure the total orbital phase to ~ 10~ 3 over the 
LIGO bandwidth [l6| . these perturbation expansions must be calculated to high order. If we consider a circular orbit 
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■ in the adiabatic approximation, the signal phase $(w) is related to the orbital energy E{uf) and the radiated power 
P{uj) through the relation d 2 Q/du> 2 ~ (dE / du) / P . An accuracy of ~ 10~ 3 in the cumulative orbital phase, over the 
LIGO bandwidth, can be achieved if the perturbation expansion is calculated to 0(v e ) beyond Newtonian dynamics 
i.e., at third post-Newtonian order (3PN) This implies that we need to know E(co) and P(lo) to at least 3PN. 

Since the conservative dynamics, described in our formalism by a Lagrangian, gives E{oj), it also must be known to 
3PN. There is also a need to compute the perturbative expansion to high order, to allow numerical studies of the 
binary inspiral to be compared with the analytic PN expansion [l8| . 

Recently, Goldberger and Rothstein introduced an effective field theory (EFT) for nonrelativistic gravitational 



systems, known as NRGR [19|, l20| (see |2l| for a pedagogical introduction). EFTs are particularly well suited to 
problems with multiple physical scales, and the binary inspiral which we consider here, is one such problem. The 
hierarchy during inspiral takes the form r c <C r <C A, where r c is the radius of the compact objects, r is the 
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orbital separation, and A is the wavelength of the emitted GWs. One can use this hierarchy to set up a tower 
of EFTs which systematically account for effects at each scale. This approach disentangles the physical effects of 
different scales, resulting in more tractable calculations. With the definite power-counting scheme established by 
Goldbcrgcr and Rothstein , the EFT treatment yields a completely systematic description of the binary inspiral 
problem. In particular, divergences which arise by using point particle sources to represent the compact objects are 
well understood in a field theory setting, and the effects due to the spatial extent of the compact objects can be 
systematically parameterized by subleading operators. 

State-of-the-art calculations using traditional methods have yielded the conservative equations of motion and the 
orbital energy of a binary with spinless constituents to 3PN order 2^, 23|, 24, 25|. The 2PN dynamics of a three-body 
system have been obtained in [261 ] . but the four-body dynamics at 2PN are not known in closed form, see, for example, 
Appendix D of 27[. Furthermore, the power emitted by a binary with spinless objects through GWs is known at the 



3.5PN level 28[ or 0(v 7 ) beyond the leading quadrupolc formula. While the EFT approach lags behind the traditional 
post-Newtonian methods in terms of high precision calculations of GW observables for binaries without spin, it has 
led to a number of interesting results. Spin was incorporated into the EFT framework and the next-to-leading order 
(NLO) spin-spin potential, which enters at 3PN, has been calculated for the first time 29J. Dissipative effects such 
as absorption by BH horizons have also been considered j^]. In [3l|], the thermodynamics of compactified black 
holes were studied, and the Einstein-Infeld-Hoffman (EIH) Lagrangian and the quadrupole formula were derived in 
arbitrary dimensions in [32I ]. The EFT formalism was extended to the case of extreme mass-ratio binaries where 
the leading order (LO) self-force equation was derived 33j. Recently, it has been used in calculations which include 
interactions beyond Einstein's general relativity 34l. l35|. 

Within the EFT method, Kol and Smolkin (KS) suggested that a temporal Kaluza-Klcin paramctrization of the 
metric improved the calculational efficiency of NRGR. The KS parametrization has been shown to reduce the 
complexity required to calculate the EIH Lagrangian 1371 . Similar simplifications were observed in the computation of 
thermodynamic properties of compactified black holes [36j , and in the calculation of the NLO order spin-spin potential 

In this paper, we report the NRGR calculation of the 2PN interaction Lagrangian for a binary system with non- 
spinning compact objects. Existing EFT calculations of binary observables have been obtained at NLO for potential 
interactions and at LO for the radiated power. This work presents the first next-to-next-to-lcading order computation 
of a GW observable with the EFT method. Since the complexity of our calculation is more involved, we clearly want 
to find an optimum method to perform high precision calculations within NRGR. Previously, it was not clear how 
useful the KS parametrization would be as a computational tool beyond NLO. We address this issue in detail, and 
we show how the KS parametrization simplifies our 2PN calculations. The methods of EFT are used to systemati- 
cally determine all relevant Feynman diagrams which contribute at 2PN. In their evaluation, we encounter Feynman 
integrals corresponding to one-loop and two-loop integrals, which are computed using standard techniques. Our work 
demonstrates how the EFT method to can be used to efficiently compute conservative dynamics at high precision in 
the PN expansion. 



II. SETUP 



Here, we outline the ingredients of our calculation of the 2PN potential and discuss the simplifications we will 
employ. In this section, we do not repeat in entirety the setup of the EFT description of the binary inspiral problem, 
but refer the reader to 19 1, [if! ], and 21 1. Since we deal with conservative dynamics, we only need potential modes, 



and can simply set the radiation modes to zero. 

A. Two-body action 

The purely gravitational action is the usual Einstein-Hilbert action 

Seh = ^ m2 pi / d i x^/~^gR 1 (1) 
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where our conventions arc m 2 Pl = 1/32-irG, t]^ = diag [1, -1, -1, -1], W va p = d a Y^p - dpT» a + T^T^ - T^r^, 
and i? M „ = R a In addition, we must also fix the gauge. Our choice is harmonic gauge 
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where T M = Y^p9 al3 ■ ^ differs from the linearized harmonic gauge condition used in [19( and yields different expressions 
at the 2PN level. The advantage of using harmonic gauge is that we can compare intermediate gauge-dependent results, 
such as the equations of motion, with the existing literature 14 1. 

For the gravitational coupling to two massive compact objects, the worldline action for our binary system is given 

by 
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Here, proper time is given by dr 2 = g^dx^dx" and it is convenient to use coordinate time t to parameterize the 
worldline. The dots denote sublcading operators encoding finite size effects. Wc will ignore finitcsizc effects since they 
first enter at 5PN [19 1. 



B. Kol-Smolkin variables 

The expansion of the metric around flat Minkowski space in the weak field limit is commonly parameterized in a 
Lorentz covariant form such as 

<Jiiv = + h^v/mpu (4) 

where is the excitation around Minkowski space. The expansion of Eq. ([4]) was used in the original setup by 
Goldbcrgcr and Rothstein [l9| . Instead, we will use an alternative parametrization based on a temporal Kaluza-Klein 
decomposition, as suggested by KS [3f|. In terms of the KS variables, the metric reads 



-e 



where 7^ = (5ij + (Tij/mpi). In four dimensions, the metric excitations are now described in terms of a scalar field 
</>, a three- vector field Ai, and a 3 x 3 symmetric tensor field er^-, all of which have been normalized to have mass 
dimension 1. In the ground state, (</>) = (Ai) = (cry) = 0, and the metric reduces to Minkowski space-time. In terms 
of the KS variables, the EH action has the simple form 

Seh = -2m 2 Pl J d 4 xV7 \-R[l] - \e^F l3 F^ + 2(d i ^)(d i ^) + ...), (6) 

where Fij = diAj — djAi and the dots denote terms with time derivatives. In Eq. ([6]) all spatial indices are lowered 
or raised with 7^ or r y VJ , respectively. Some of the terms involving time derivatives are needed for our calculation 
along with the expansion of the gauge fixing action of Eq. ([2]), and their contributions enter the Feynman rules given 
in Sec. EH 

Given the gravitational action Eqs. |T]) and j2]), the worldline action from Eq. ((3|), and the KS metric parametriza- 
tion from Eq. |5]), we can proceed to derive Feynman rules for the gravitational self-couplings and the couplings to the 
compact object worldlines. The procedure is equivalent to the one in [l9| . but is now implemented in terms of the KS 
fields. For the gravitational self-couplings and worldline couplings, there are infinitely many terms in the respective 
actions. To avoid unnecessary calculation at a given order, a method of systematic computation must be established. 
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FIG. 1: Order G topology. The single solid line denotes a generic graviton field, either <j>, Ai, or <Tij, and double solid lines 
denote the worldlines of the binary constituents. 




(a) (b) 
FIG. 2: Topologies at order G 2 . 



C. Power counting and Feynman diagram topologies 

The power-counting rules for our calculation arc the same as those in the original EFT setup developed in [l9| . 
However, since there are no radiation fields in our 2PN potential calculation, we use a simplified procedure where we 
power count our Feynman diagrams relative to the LO Newtonian potential. For a bound state, the virial theorem 
relates the orbital velocity v and Newton's constant G via, v 2 ~ Gm/r, and as shown in [l9| . one can power count 
all contributions in powers of the orbital velocity v. Keeping the virial theorem in mind, we count powers of G and 
powers of v 2 separately. For the power counting we do not keep track of factors of mass m and separation r associated 
to each G. These will be generated appropriately in the calculation of the diagrams. In this scheme, the diagram 
responsible for the Newtonian potential scales as O(Gv ). For the 2PN potential calculation, we need to include all 
diagrams which scale as 0{Gv A ), 0(G 2 v 2 ), and O(G 3 v ). 

To determine the relevant Feynman diagrams at 2PN, the first stage is to generate all relevant diagram topologies. To 
do this, one first counts in powers of G, since all interaction vertices will scale with a power of G (recall mpi ~ G -1 / 2 ). 
Powers of velocity are inserted later, when all relevant diagram topologies have been established. There are two rules 
to consider when counting a topology in factors of G. Firstly, when there are n gravitons attached to a worldline, this 
component receives a factor of G n l 2 . Secondly, each n-graviton self-interaction vertex receives a factor of G™/ 2-1 . 
Note that a propagator does not receive any factors. Using these rules any diagram topology can be counted in powers 
of G. 

If we now proceed to apply these rules, we can only have a single topo logy at 0(G), as shown in Fig. [T] All 
topologies which simply renormalize the mass are omitted, as discussed in [19(. Diagrams which involve a graviton 
line starting and ending on the same worldline, without any intermediate interaction, fall into this category. At the 
next order, 0(G 2 ), we have two topologies, as shown in Fig. [2] Since the worldlines are static sources and do not 
propagate, the diagrams in Fig. [2] do not involve any loops. Diagrams with gravitational loops yield quantum effects 
and are therefore ignored, and since only massless fields propagate in NRGR, the expansion in loops corresponds to 
an expansion in powers of h [3^|. At 0(G 3 ), the total number of possible topologies is nine, as shown in Fig. [3] These 
diagrams are the relevant topologies for the static component of the 2PN potential. In Figs. [2] and [3l we have not 
drawn the topologies with the two worldlines interchanged, although they will be required for our calculations. 

Now that we have established all topologies relevant to the calculation at 2PN, the powers of the orbital velocity v 
must be counted. There are two possible sources for factors of v in our calculations: (1) fields coupling to the worldlines, 
where the LO couplings of <fi, Ai, and cr^ are (D(v ), Oiv 1 ), and 0(v 2 ), respectively, and (2) time derivatives, where 
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FIG. 3: Topologies at order G 



3q ^ v/r for potential modes. Time derivatives can arise from either purely gravitational interaction vertices or from 
propagator insertions, where each propagator insertion counts as 0(v 2 ). 

With all relevant topologies and the counting rules for v, the final diagrams at 2PN can now be constructed. Recall 
that the 2PN diagrams scale as either 0(Gv 4 ), 0(G 2 v 2 ), or 0(G 3 v°). The diagrams are generated by first populating 
the topologies with the three gravitational fields (f>, Ai, and Uij in all possible combinations. We then count the factors 
of v which result from the gravitational fields coupling to the worldlinc, and any time derivatives acting on internal 
vertices. Propagator insertions are also counted where appropriate. The powers of v from the worldline vertices are 
kept to all orders in the calculations and only expanded to 2PN in the final result. The diagrams relevant at 2PN 
which result from this procedure are shown in Sec. Figs. IU03 and[6l 



D. Advantages of the Kol-Smolkin variables 

In the KS parametrization of Eq. ([5]), we have introduced three fields, <fr, Ai, and er^-, instead of the standard h^ v 
paramctrization of Eq. ([4]). It will be necessary to keep track of these new fields, so one should ask what has been 
gained from a calculational perspective by using the KS parametrization. This question is answered by considering 
the static diagrams, i.e. those with no velocity factors. First, consider the topology Fig. &b) which involves a three- 
graviton vertex. This was used in [l^ to construct the EIH potential in the variables. When the KS variables 
are employed, it was shown that this topology docs not contribute at 1PN order We discuss this in more detail 
since it points us to the source of the advantages of the KS variables. 

When working in the parametrization, one often considers its components Hqq, hoi, and hij separately, since 
their LO coupling to the worldline scale as O(v ), OfV), and 0(v 2 ), respectively. So in practice, in the /i M „ and KS 
parametrizations, one must keep track of the same number of components or fields. Returning to our discussion, if 
we consider static diagrams, only the hoo component can couple to the worldlinc. Let us consider the static limit of 
the first topology which arises with gravitational self-couplings, Fig. [2^6). With hoo coupling to each worldline at 
LO, without any propagator insertions, the only possible source of powers of v would be the three-graviton vertex. 
Naively, one would isolate the action component with three powers of hoo, which yields the 3-/ioo vertex. Interestingly, 
one finds that the corresponding term in the action is proportional to hoo(dohoo) 2 , meaning that the 3-hoo vertex has 
a power of v 2 associated with it. So does this mean that the topology in Fig. \2lb) does not enter at 1PN in the h^ v 
parametrization? If this were true, then the EIH calculation of [19j would be contradicted. The reason why the above 
conclusion is wrong is the presence of mixing of the h^ v components. 
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To explain this point, consider the h^ u graviton propagator in harmonic gauge, which is given by (Th^hap) ~ 
(j)ticJ]v(3 + f]ni3f]va — VftvVap)- Examining the tensor structure of the propagator, one can see hoo mixes with the trace 
of hij, since (Thoohij) ~ 5ij. Thus, a graviton which starts as hoo at the worldline coupling can be either hoo or ha 
at the three-graviton vertex. So instead of the 3-/ioo vertex, one has to use the three-point function (T/ioo^-oo^oo)- 
This three-point function is obtained by contracting the 3-h^ vertex with three graviton propagators and setting all 
remaining free indices to 0. One then finds that there are components of the three-point function (T/ioo/ioo^oo) which 
do not involve time derivatives, and so the topology Fig. \^b) is required at 1PN in the h^ v parametrization. Clearly, 
the topology Fig. ^b) only survives at 1PN in the h^ u parametrization because h 00 can mix with ha as it propagates. 

The major advantage of the KS parametrization in harmonic gauge is that it removes this kind of mixing, since the 
two-point functions between the three fields are zero: (Tcf>Ai) = (T<f><jjk) = (TAiUjk) = 0. For the topology Fig. 
to contribute to the EIH potential in the KS parametrization, all of the fields coupling to the worldlines must be tfi's, 
since only they couple at O(v ). Knowing that <f> cannot mix with Ai or cry, the 3-(f> vertex must be used. When the 
gravitational Lagrangian is examined, the only relevant term is proportional to 0</> 2 , where the two time derivatives 
in this term cause the 3-<f> vertex to be of 0(v 2 ) in the orbital velocity. Therefore, the diagrams resulting from the 
topology Fig. [21J6) do not contribute until 2PN in the KS parametrization. 

This argument can be extended to the n-<f> vertex (n > 2), where the relevant term in the gravitational action is 
~ exp(— A(j)/mpi)(j) 2 . Any n-<fi vertex is then suppressed by one order in the PN expansion because there are always 
two time derivatives. Thus, in the KS parametrization, the 0(G 3 ) topologies Fig. (c), (/), and (g) will not enter 
at 2PN. This will reduce the calculational effort considerably Also, the fact that the topology Fig. first enters 
at 2PN, rather than 1PN, means that propagator insertions are not required for this topology at 2PN. 

We have established that the KS variables in harmonic gauge provide important calculational advantages over 
the h^y parametrization used in [l^| at 2PN order. At 1PN, this advantage was modest, since only one diagram 
topology was eliminated. At 2PN, we have found that four topologies at O(G 3 v ) were removed, and one did not 
need propagator insertions at 0(G 2 v 2 ) on the topology Fig. &b). Clearly as we proceed to higher order in the PN 
expansion, the advantages of the KS variables become more important. 

III. CALCULATION OF FEYNMAN DIAGRAMS 

In this section, we present the calculation of the Feynman diagrams which are required to determine the 2PN 
interaction potential between the binary constituents. Using the Feynman rules presented below, each diagram equals 
—i J dtV from which we extract its contribution to the potential V. The symmetry factor for each diagram is 
computed in the usual manner [4oj |. but one needs to account for the nonpropagating classical sources. As discussed 
in Sec. Ill CI there are three classes of diagrams with the following power counting: O(G 3 v ), 0(G 2 v 2 ), and 0(Gv 4 ). 
We discuss how to evaluate each in turn, starting with the simplest at C(Gv 4 ), and finishing with the more complex 
diagrams at O(G 3 v ). 

When calculating the diagrams, we work with the spatial Fourier transform for the gravitational potential modes. 
For example, for the </> field, we will work in terms of 4>k(t), which is the spatial Fourier transform of the coordinate 
space field 4>(t, x) = J k e lk ' x </>k(i), where J k = J . We also make clear that after expanding the metric in the KS 
fields <j), Ai, and cry, we do not distinguish between upper and lower spatial indices; they are lowered and raised with 
5ij and 5 V , respectively Lastly, since we are interested in the long range potential only, we drop irrelevant contact 
terms, such as V ~ S(r), wherever they appear in the calculations. 

A. Order Gv 4 diagrams 

As shown in Fig. [4j there are six diagrams to be evaluated at 0(Gv 4 ). The contributions from the diagrams (a), 
(d), and (/) are the easiest of all diagrams to compute at 2PN, since we only need to compute simple one-graviton 
exchange diagrams. 

The Feynman rules for couplings to the worldlines are derived from the point particle action in Eq. ([3]), which is 
expanded in the gravitational fields to the required order for a given Feynman rule. The Feynman rules for a <ft, Ai, 
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FIG. 4: Order Gu 4 dia grams at 2PN. Here the dashed, wavy, and double wavy lines, represent the <j>, Ai, and cnj fields, 
respectively. A cross denotes a propagator insertion. 



and <Tj,- coupled to the worldline are 



mpi 



dt I e 

k 



ik-x 



VT- 



(7) 



= inL[ dt [ e i** 

mpi J A 



|k^ = / dt I <' kx V ' Vf . (9) 



2topi V A \/l 



-v 



The double solid lines are the worldlines, and the dashed, wavy, and double wavy lines represent the <p, Ai, and tr^ 
fields, respectively. Note that these rules are exact in v. We will expand in v at the end of the calculation to isolate 
the required terms at 2PN. 

To compute diagrams (a), (d), and (/), we now need the propagators for each gravitational field. From the purely 
gravitational action, the sum of Eqs. |T]) and @, the propagators can be derived. These potential propagators are 
nonrclativistic and instantaneous, and for the </>, Aj, and gravitational fields we have, respectively, 

(7>p(i Q )(Mi b )) = -\(^) 3 S 3 (p + q)A I 6(t a - t b ) (10) 

(TAi(t a )A{(t b )) = i(2^) 3 <5 3 (p + qf-^Sita t b ) (11) 



(Ta»(t a )a%(t b )) = -(27rfS 3 (p + q ) l -^-S(t a - t b ), (12) 



where P^' M = ±(5 ik 5 jl + S a 6^ k - 2S ll S kl ). 

The three single-graviton exchange diagrams are easily computed with the above propagators and worldline cou- 
plings. At 0(Gv 4 ), the only integral in k which arises is a Fourier transform, which is evaluated using the d-dimensional 
formula 

d^ 1 1 T(d/2-a) /-2N"-"/ a 



(27r) d (k 2 )" (47r) d / 2 T{a) 



^ = ' hr ■ (13) 
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Doing the Fourier transform, and extracting the 2PN component, gives the following potential contributions: 
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where the notation we use here and subsequently is r = xi(f) — X2(i), r = |r|, and n = r/r. The labels 1 and 2 
are used for the left and right worldlincs, respectively. Although r, Vi, etc., depend on t, we have suppressed this 
dependence above. We will continue to do this when convenient. 

Diagrams (b), (c), and (e) all involve propagator insertions. These insertions account for the corrections to the 
instantaneous nature of the nonrelativistic potential propagators in Eqs. (fT0|) - (|T2|) . These are included systematically 
as a perturbation, where each propagator insertion is suppressed with respect to the propagator by a power of v 2 . 
The insertions are obtained from the terms in the action which are quadratic in the fields and which involve time 
derivatives. The required propagator insertions at 2PN are derived to be 



- -X - - 



^(27r) 3 6 3 (p + q)^d ta d tb S(t a -t b ) 
I(27rf6 3 (p + q)^dldl6(t a -t b ) 



i 

P 

2tt^ 3 
2 K ■ p 



= -(27r) 3 5 3 (p + q)— d ta d tb S(t a - t b ), 



(17) 
(18) 
(19) 



where one cross denotes a single propagator insertion and two crosses denote two propagator insertions. Here, we have 
chosen the simplest symmetric form of the time derivatives; however, other choices are possible. Note that any other 
choice of time derivatives yields potentials which can be related to potentials computed with the simplest symmetric 
choice by a total time derivative, and are thus physically equivalent. 

To demonstrate how the diagrams at 0(Gv 4 ) arc computed, we show how to calculate diagram (b) in some detail. 
The first step is to form the expression of the diagram. This is achieved by coupling Eq. (fT7|) to two copies of Eq. ([7|) 
for worldlines 1 and 2. The diagram is given as follows, where we have done the q integration to eliminate <5 3 (p + q), 



dtV 



8m 2 Pl 



dt n dt h 



e ip-( Xl (t a )- X2 (t b )) (1 + y 1 (f a )2)( 1 + y 2 (t b ) 2 ) 

P 4 (l-V 1 (t a ) 2 ) 1 /2(l_ V2 ( 4 )2)l/2 



dt a d tb S(t a -t b ). 



(20) 



Care must be taken when labeling the time for each worldline Feynman rule, and here we have used t a and t b . 
Since each worldline interaction is a different Feynman rule, it must therefore have a different dummy variable. 
Before explicitly using the time derivatives, we employ the Fourier integral from Eq. (|13p to compute the remaining 
momentum integral. Now, an integration by parts on both time derivatives acting on S(t a — t b ) is required, which 
gives, 



dtV = 



—iGm,\m,2 



dt a dt b 5(t a - t b ) \ d ta d t 



|xi(t a ) - x 2 (t b )| 



(l+v 1 (^) 2 )(l+v 2 (t b ) 2 ) 

(l-V 1 (t a ) 2 )V2(l_ V2 ( 4 )2)l/2j 



(21) 



At this stage, the time derivatives must be allowed to act on the square bracket. Only after this point can the 5(t a —t b ) 
be used. It is clear that acceleration-dependent terms will arise. The potential contribution we obtain is exact to all 
orders in the PN expansion. Since we are only interested in the 2PN potential, we extract this piece from the result. 
For this purpose, we note that accelerations are power-counted as o ~ v 2 /r. The final contribution, V^ b \ is given 



9 



(fj) 



(h) 



U) 



(k) 



-< 



— v 



--V 



(0 



(m) 



n 



(p) 



FIG. 5: Order G 2 v 2 diagrams at 2PN. 



below, along with the other one-graviton exchange diagrams with propagator insertions at 0(Gv ) 
W = - (f ( Vl • v 2 - n • Vl „ • v 2 ) (v? + v- 

- Gmim 2 Qa x ■ ■ v 2 - ^a 2 ■ v 2 n ■ ViJ (22) 
V^ c) = ( g v I v 2 + ^(vi • v 2 ) 2 + -(n • Vl ) 2 (n • v 2 ) 2 

- -Vi • v 2 n • vin • v 2 - ^v 2 (n • v 2 ) 2 - ^v 2 (n • v : ) 2 
z o o 

- Gmim 2 Qai • n (v 2 - (n • v 2 ) 2 ) - ^a 2 • n (v 2 - (n • Vl ) 2 ) 

1 1 
+ -a : • v 2 n • v 2 - -a 2 • • Vi 

- Gmxm^r ( — -ai ■ a 2 — -ai • na 2 ■ n ] (23) 



- Gmim 2 (— 2ai • v 2 n • v 2 + 2a 2 • vin • Vi) 

- Gm\m<zr (2ai • a 2 ) . (24) 

Note that all the potential contributions V^ b \ V^ c \ and V^ e \ contain acceleration-dependent terms. They are 
generated from either a time derivative acting on a worldline velocity factor or from two time derivatives acting on 
xi or x 2 . 

B. Order G 2 v 2 diagrams 



The first three diagrams at 0(G 2 v 2 ), (g), (h), and (i) of Fig. [5j are obtained by modifying the 0(G) topology 
through the addition of one extra gravitational leg connecting the two worldlines. As we will see, such diagrams, 
where a topology of lower order in G is augmented by an additional leg between the worldlines, are particularly 
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simple to compute because their expressions factorize. To compute these diagrams we need two additional Fcynman 
rules for the worldline couplings. The first has two <f> fields coupling to the worldline, and the second has a 4> and an 
Ai coupling to the worldline. These are respectively, 



m 



PI J A,q (1-V 2 ) 



'PI 



(It 



^(k+q)-x ( 1 - 3v2 ) V » 



k.q 



(1-v 2 ) 



3/2 



(26) 



These rules are again exact in the orbital velocity v. 

We demonstrate how to calculate the three diagrams (g), (h), and (i), by considering diagram (h). One can write 
this diagram by using Eqs. ||7J), (JSJ), and |26|) for the worldline couplings, and Eqs. (fTO)) and (fTTj) for the propagators. 
The diagram is 



m ^ m 2 I M _M ldt XU. . _*.W* . / p<[(k+P)-xi(t.)+l.x a (t i )+q.x a (t.)] 



iity= ' 4 / dt a dt b dt c 5{t a - t b )S(t a - t c ) / e 
lbm p; J Jk,i,p,q 

(2^) 3 £ 3 (k + 1) (2tt) 3 £ 3 (p + q) vifa) ■ v 2 (t c )(l - 3 Vl (t a ) 2 )(l + v 2 fa) 2 ) 
X k2p2 (l-v 1 (t a ) 2 ) 3 / 2 (l-v 2 (4) 2 ) 1 / 2 (l-v 2 (t c ) 2 ) 1 /2 



(27) 



The symmetry factor of this diagram is 1. Since there are no time derivatives in this expression, the time delta 
functions can be used immediately. We can also use the three-momentum delta functions; doing so leads to 

rniimj f f e* r Vl ■ v 2 (l - 3v 2 )(l + v 2 ) 
" WmU J d Vk. P k 2 p 2 (1-^(1-^) ■ ^ j 

Looking at this equation, it is clear that the two momentum integrations factorize into two Fourier transforms. They 
are seen to be the same in this diagram, and we can simply use the Fourier integral, Eq. (|13[) . Our final result for 
the exact potential contribution from diagram (h) to all orders in the PN expansion, is then 



r (h) _ -4G 2 m 1 m 2 v 1 ■ v 2 (l - 3v 2 )(l + vj ) 

(1-V 2 )3/2(1~V 2 ) 



K Exact — „2 ..2\S/2/'i ,.2\ 1 V. 



The calculation of the other two diagrams with this topology, diagrams (g) and (i), proceed similarly. However, 
when propagator insertions are present, as in diagram (i), the time delta functions cannot be used immediately. One 
must first integrate by parts, and act with the time derivatives, before using the time delta functions. The procedure 
is analogous to diagram (b) in Sec. IIII AI At 2PN, we have the following contributions to the potential: 

V«-^(-|vf + |v5) (30) 
^)=^=l ( -4 Vl .v 2 ) (31) 

j^v4 n v2)2 + i ^ ^ _ n \ (32) 

Although we do not display the diagrams where the two worldlines arc interchanged, they must be accounted for in 
the final 2PN Lagrangian. 

We now proceed to the diagrams where the topology is that of Fig. [DJb), where a three-graviton vertex mediates 
the interaction. First, the vertices for the gravitational self-interactions must be derived, and this is done in a 
straightforward way from the sum of the Einstein- Hilbert action, Eq. |T]), and the gauge fixing action, Eq. At 
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2PN, the KS parametrization simplifies the vertex derivations, since we need at most one leg to be the tensor field 
Uij. Once the appropriate terms in the action are established via power-counting factors of v, the Feynman rule for 
the vertex is computed. Then each leg of a given vertex is multiplied by the appropriate propagator to obtain the 
three-point function. For diagrams (J), (k), (/), (to), (n), (o), and (p) of Fig. we use the following three-point 
functions: 

{T<f> a( j> b( f> c ) = Yg^T) ( dt « 6 ( ta ~ *°)) - *c)) + (dtj(t a - 4)) (d t J(t c - t b ) 



d tb S(tb - t a )) (d t J(t c - t a )) J> x j-^p-j (2tt) 3 6 3 {k a + k b + kc) (33) 



(Td> a <t> b Ai) = - jg^-| (kjft. + Kdt b )s(t a - 4)5(4 - t ) 1 x (2^) 3 <5 3 (k a + k fc + kc) (34) 




k^k? — k'k* + ki • k fiv 

(T^ a AlAi) = — S(t. a - t c )6(t b - t c ) x -*-£ ^ 2 ^ 2 c (2^) 3 <5 3 (k a + k 6 + k c ) (35) 

(Tfafo***) = 8(t a - t c )5{t b - tc) x !^-|+MM ( 27r ) 3 <5 3 (k Q + k fa + ke), (36) 

loTOpz K a K^K c 

where the lower index on a field denotes its dependence on t and k. 

Because of the nonlinear structure of the diagrams (j) through (p), their contributions do not factorize, and we 
have to calculate an integral which corresponds to a one-loop Feynman integral. To demonstrate the calculation of 
this class of diagrams, we consider diagram (fc). This diagram uses the three-point function (T(f) a A z b A J c \ from Eq. 
([55)) . The diagram is formed by coupling Eq. (|3"5|) to the worldlines using Eqs. ([7]) and © and the symmetry factor 
of this diagram is 1. Given there are no time derivatives, we are free to do two of the time integrations immediately. 
We also do the integration over the three-momentum which labels the <f> leg. The result of these operations is 

MV 4m 4 p; J ^(l-vf)i/2(l_ v 2) 7 W 
where the integral I lJ (r) is given by 

7p, q p 2 q 2 (p + q) 2 i P p 2 7 q q 2 (q + p) 2 

The next step is to evaluate the momentum integral, /^ (r). From Eq. (f3"5| it is apparent that 7 y (r) is calculated by 
first computing the q-integration which corresponds to a one-loop Feynman integral and subsequently performing a 
Fourier integration over p. The Feynman integral is of vector nature, and we reduce it to a scalar integral in the usual 
way by noting that its expression must be proportional to p a . The scalar integral is then computed in dimensional 
regularization using the d-dimensional master integral 

d d k 1 1 T{ni+n 2 -d/2)T{d/2-n x )T{d/2-n 2 ) 2dn _ ni _ n ^ 



(2n) d [(k + p) 2 ]" 1 [k 2 ]" 2 (4^) d / 2 r(m)r(n 2 ) T(d-nx-n 2 ) 
The resulting expression for the integral (r) then becomes 



S lj [ ■ 1 

r,(r) = ~i6 I etp ' r —> (40) 



p 



and the remaining Fourier transform integral is performed with Eq. (| 13[) . Putting all terms back together gives the 
exact potential for diagram (k), 



r(fe ) _ G 2 m 1 ml 8vi-v 2 (l + V2 

(1-V 2 )V2(1_ V 2)' 



v Exact — <i fl ..2M/2C1 V^V 



The 2PN piece of this potential is extracted below. 
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FIG. 6: Order G 3 w° diagrams at 2PN. 



The remaining diagrams from (j) through (p) are all evaluated using the same methods as explained for diagram (fc). 
The only complication are the time derivatives which occur in the three-point functions (T</> a 0b</> c ) and (T<j> a (l>bA l c }, 
Eqs. (|33[) and (J33J), respectively. So when calculating the diagrams (j), (m), and (n), it is necessary to integrate by 
parts as explained in Sec. IIII Al We use Eq. (and its vector or tensor extensions) for any one-loop integrals, 

and Eq. <| 1 3j) for the Fourier transforms. Additionally, for some diagrams we need vector and tensor Fourier integrals 
which are obtained from the scalar Fourier integral in Eq. (fl~3"|) by taking derivatives with respect to r. At 0(G 2 v 2 ), 
the contributions from the individual diagrams with the topology Fig. [2fb) are 

(42) 

(43) 
(44) 
(45) 
(46) 
(47) 
(48) 



yO') 


G 2 mim\ j 

r 2 \ 


jv2-i(n.v 2 ) 2 


— 2vx • v 2 + 4n • vxn • v 2 


v (k) 


G 2 m\m\ 
r i \ 


8vi ■ V2) 




yW 


G 2 mira2 , 

r 2 V 


-4v 2 2 ) 




y(m) 


G 2 m\m\ ( 


2w\ -4(n- v 2 ) 2 + 


2vi • v 2 — 4n • vin • v 2 ) 


y(") 


G 2 raim 2 , 
r 2 \ 


— vx • v 2 + n • vin 


• v 2 ) 


y(°) 


G 2 mim 2 ( 


-2 V 2 + 4(n-v 2 ) 2 ) 




yip) 


G 2 niim 2 i 

r 2 \ 




)■ 



All of the above contributions with (1 <-> 2) must be added to the final potential in order to account for the diagrams 
with the two worldlincs interchanged. 



C. Order G 3 v° diagrams 

Finally, we progress onto the diagrams at O(G 3 v ). The diagrams which must be evaluated are (q), (r), (s), (t), 
and (u), as shown in Fig. [SJ For diagram (q) a new worldlinc vertex will be required coupling three </> fields to the 
worldline at one point. The Feynman rule is 

,( k+q+ pVx l + llv 2 + Hv 4 +v 6 f4Q) 
, P (1-v 2 ) 

which is exact to all orders in v. At this order, the potential contribution at 2PN is static and there are no velocity 
factors in the final potential contributions. Given this is the case, we do not keep any velocity factors from the 
worldlinc vertices when we discuss the calculations below. 

The two simple diagrams at this order are (q) and (s). These diagrams are easy because there are no internal 
vertices - only propagators and worldlinc couplings. The diagrams therefore factor into simple Fourier transforms, 
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which are evaluated using Eq. (| 13|) . This is analogous to the evaluation of diagram (h) in Sec. IIII Bl The symmetry 
factors are 1/6 and 1 for (q) and (s), respectively. Computing each diagram gives the following potential contributions: 

v(q) = GhnM 

^=- G y , (si) 

where the contribution of diagram (q) must be added with (1 <-> 2) in the final potential. 

The next set of diagrams at 0(G 3 v a ) are no longer as simple as (q) and (s). This is due to the presence of the 
four-point function, which is used to construct the diagrams (r), (t), and (u). The only four-point function needed 
at 2PN is derived by joining two (fxpcjij vertices with an intermediate propagator. Eq. (|12p . and attaching <f> 
propagators, given in Eq. (flO|) , to the external legs. This is given by 



T(f>a4>b §c4>d ) = — 5- ${ta ~ t d )5(t b - td)8(t c - t d ) 

k a • k c k ■ k^ + k a ■ k^kb ■ k c — k a • k k c ■ 3 3 /m\ 

X (k Q+ k fc)2 kgkg k 2 k 2 (2-) * + k h + k c + M, (52) 

where each set of two (f> fields connected by a brace have the same intermediate (jxf)<7ij vertex. 

Of the three diagrams left to compute, (t) is the most involved, so we will compute it explicitly. We first note 
that the symmetry factor of this diagram is 1/2. The diagram is constructed in the usual manner, using Eq. ([7]) 
and Eq. (|52p . but we will only work to 0(G 3 v a ) in what follows neglecting all velocity factors of the vertices. First, 
we integrate over all the time delta functions. Then, we relabel the three- momentum of the four-point function as 
(k a , k(,, k c , krf) — > (k, q, p, 1) and couple k and p to worldline 1 and q and 1 to worldline 2. The 1 integration is then 
performed using the (5-function, which takes 1 — > — (k + p + q) . After some algebra, redefining of k — > k + p is seen 
to be useful. This gives the following expression: 

dtv= *mfm| rr elk . r + + + 

256m° pl J J k 

where the terms ^i(k), In(k), Jni(k), and iiv(k) are the integrals 

7 x (k) = f -k.pq 2 /Q(k,p,q) (54) 



7 n (k) = / k.q P 7Q(k,p,q) (55) 
-'p,q 

7„l(k) = [ p 2 q 2 /g(k,p,q) (56) 
-'p.q 

7 IV (k) = j -k 2 p-q/Q(k,p,q), (57) 

and Q(k, p, q) = (k + q) 2 (k — p) 2 (k + q p) 2 q 2 p 2 . These integrals correspond to two-loop Feynman integrals, and 
we will show how to compute them efficiently by reducing them to the one- loop master integral in Eq. (1391) . 

The first three integrals ii(k), iii(k), and iin(k) can be computed straightforwardly using the d-dimensional one- 
loop scalar integral from Eq. (|39|) twice, on both the p and q integrations. One can see that 2i(k) = /n(k) by 
making the replacement (p,q) — > (— q, — p), in iii(k). The reason why one can evaluate these two- loop integrals 
employing solely a single one-loop master integral, is that these integrals correspond to nested loop diagrams. Here, 
it is important to keep the dimension d of the integrals arbitrary and only set d = 3 after the final Fourier integration 
over k has been performed. 



The remaining two- loop integral, iiy(k), requires an integration by parts trick 4l| in order to reduce it in terms of 
our master integral of Eq. (|3"9"|) . But first, we redefine q — > q + k. This allows us to split the integral into four terms, 
three of which can be evaluated using same methods as discussed in the previous paragraph. The remaining term is 
a two-loop integral with five factors in the denominator 

k 4 f 1 



2 7 p , q q^ (q - k)2 (q + p)2 (p + k )2p2 ' ( 58 ) 
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Using the integration by parts trick we can write this integral as 

1 2 



q 2 (q-k) 2 (q + p) 2 (p + k) 2 p 2 d - 4 L_ VpV(q + k) 2 (p + k)* p 2 q 2 (p + q) 2 (p + k)* 



(59) 



Once in this form, we can evaluate it using Eq. (|39p twice. 

After all the Feynman integrals have been evaluated, we perform the remaining Fourier transformation integral 
using Eq. (|13p . Subsequently, we can set d = 3 and find a contribution to the potential 

v{t)= _2G*rn\ml (60) 

The remaining two diagrams, (r) and (w), are simpler to compute. They are evaluated with the methods we have 
discussed, where once again it is necessary to compute two-loop integrals by applying the master integral Eq. (|39p to 
two integrations. The results of the calculations are 

V ir) = gnwgj (61) 

V {u) = 0, (62) 

and the contribution of the potential must be added with (1 «-> 2) in the final potential. It is worth noting that 
y( u ) = o because the diagram gives a purely short distance contribution V ~ S(r), which we have dropped. 



IV. RESULTS 



Having computed all the diagrams, we can now construct the 2PN interaction Lagrangian. There are two contribu- 
tions which we need. The first is the kinetic energy, which is obtained by expanding the matter coupling action, Eq. 
©, to 0(v 6 ) while setting all fields to zero. The final potential comes by summing each potential contribution from 
the diagrams, 

y(a) t0 y(u)_ At thig gt 

we add all contributions with the worldlines interchanged, as appropriate. 
Our 2PN interaction Lagrangian for a binary system is then given by 



L 



2PN 
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Gmi77i 2 M 4 5 2 3 2 3 22 1 2 

+ I gVj - -v x vi • v 2 - - Vl n • Vl n ■ v 2 + ^ v i v 2 + g( v i ' v 2) 

13 3 ^ 

- - v 2 (n ■ v 2 ) 2 + -n • • v 2 v x ■ v 2 + — (n • vx) 2 (n ■ v 2 ) 2 

( 1 2 3 7 1 , -.2 

+Gm 1 m 2 I -ai • nv 2 + -a x ■ v : n ■ v 2 - -ai • v 2 n • v 2 - -ai • n(n ■ v 2 ) 

(™ 1 ^ 

+Gm 1 m 2 r I — ai • a 2 - — ai • na 2 • n I 



G 2 mimi ( 7 9 „ ? 7 1, 
+ r 2 [I"! + 2v 2 - 2 Vl ' V2 + 2 (n ' Vl - 

G 3 mimn 3G 3 m?mo ,„ „. 

J 2 +(1^2), (63) 

where (1 <-» 2) refers to all terms given previously, with the labels 1 and 2 interchanged. Note that n — > n under 
this exchange. 

Upon comparison to the 2PN piece of the Lagrangian in Eq. (174) in it is apparent that our Lagrangian 
does not have the same form. However, when the equations of motion (EOM) are computed we recover the same 



2PN EOM as in 14| . This shows the physical equivalence of our Lagrangian and the one in 14[ . At the level of the 
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Lagrangian, this can be shown directly by relating our result in harmonic gauge to the standard 2PN Lagrangian in 



harmonic gauge of 1J] via a total derivative and two double zero terms 42j. A double zero term is an expression 
which vanishes at 2PN by the use of the lower order EOM. Adding such a term modifies the form of the Lagrangian. 
but does not change the EOM or the choice of gauge. The terms we add are given by 



6L1 = 



1 Gm\mi 
8 r 
15 



r • ai 



Gran 



r • a 2 



8L-2 = - —Gmim 2 r ( a* 



Gm 2 i 
— 3— r 



a 9 =— v 



6L* = 



dt 



7 G 2 mim2 



(m2r • V2 — rriiT ■ vi) 



3 Gm,\mi ii o \ 

(r • viv 2 - r • vavJ 

4 r \ i 



(64) 
(65) 
(66) 



where the first two terms are double zero terms and the last one is a total time derivative. By adding these terms 
to our Lagrangian at 2PN, Eq. (f6"3")) . we obtain the Lagrangian in 14[. We again emphasize that applying these 
transformations leaves the EOM unchanged, so that we remain in harmonic gauge. 

We have also calculated the observable E(u) for a circular orbit and recover the known 2PN result [l4(. We 
note that if the accelerations were replaced in our Lagrangian, using the LO EOM, E(uj) is unchanged, although all 
intermediate expressions are different. This is because E(cu) is physical and therefore gauge invariant, and using the 
EOM at the level of the Lagrangian corresponds to a change in gauge. 



CONCLUSION 



We have demonstrated how to use the EFT method 19( to efficiently calculate the next-to-next-to-leading order 



Lagrangian describing the conservative dynamics of a binary system. We have shown how to systematically determine 
all Fcynman diagrams which contribute to the 2PN Lagrangian through the use of power counting in G and v 2 . The 
calculation involved 21 distinct diagrams, and we encountered integrals which correspond to one-loop and two-loop 
Fcynman integrals. All diagrams could be computed with only two master integrals. This demonstrates the efficiency 
of the EFT method for the calculation of the 2PN Lagrangian. 

Instead of a usual Lorentz covariant metric paramctrization, wc employed a temporal Kaluza-Klcin parametrization 
of the metric, proposed by Kol and Smolkin 36|, 37 1 . When compared to the h^ v parametrization used in [l9| , we found 



that the KS variables reduced the number of diagrams by four at 0(G 3 v°) and avoided any propagator insertions in 
the three-graviton diagrams at 0(G 2 v 2 ). This reduced the amount of calculation significantly at 2PN. More generally, 
we conclude that the KS variables will significantly improve any calculations which arc performed at higher order 
using NRGR methods. This is due to (1) the simple propagators obtained with this gauge and parametrization choice, 

(2) the suppression of <f> n vertices by one order in the PN expansion due to two time derivatives acting on the vertex, 

(3) not requiring time insertions in the 4> n topologies at the next order in the PN expansion, and (4) the advantageous 
compounding effect of the former points when higher order expansions are calculated. 

An interesting extension of the current calculation is to compute the 3PN potential. Analogously to 2PN, where 
the four-point function formed diagrams which were the most difficult to evaluate, the hardest part of the calculation 
at 3PN is expected to involve a gravitational four-point function in O(G i v ) diagrams. Also, the number of diagrams 
one must evaluate at 3PN is much greater. This means more terms in the action will have to be computed to account 
for new vertices at O(G 4 v ), 0(G 3 v 2 ), 0(G 2 v 4 ), and 0(Gv 6 ). Our preliminary estimates show there are more than 
100 diagrams at 3PN, compared to the 21 at 2PN. At 3PN, it will again be the case that a reduction in the number 
of diagrams can be achieved if the KS variables are employed. So while the 3PN calculation may prove challenging, 
the relative ease of the calculation at 2PN indicates that the NRGR calculation of the 3PN potential can be done. 
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